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CH ABSTRACT 

We present a deep [Oil] emission line survey of faint galaxies (22.5 < Kab < 24) in the 
i Chandra Deep Field South and the FIRES field. With these data we measure the star 

formation rate (SFR) in galaxies in the stellar mass range 8.85 < log (M*/M ) < 9.5 
at 0.62 < z < 0.885, to a limit of SFR~ O.lM ?/r _1 . The presence of a massive 
cluster (MS 1054-03) in the FIRES field, and of significant large scale structure in the 
CDFS field, allows us to study the environmental dependence of SFRs amongst this 
population of low-mass galaxies. Comparing our results with more massive galaxies 
at this epoch, with our previous survey (ROLES) at the higher redshift z ~ 1, and 
with SDSS Stripe 82 data, we find no significant evolution of the stellar mass function 
of star-forming galaxies between z = and z ~ 1, and no evidence that its shape 
depends on environment. The correlation between specific star formation rate (sSFR) 
and stellar mass at z ~ 0.75 has a power-law slope of f3 ~ —0.2, with evidence for 
a steeper relation at the lowest masses. The normalization of this correlation lies as 
expected between that corresponding to z ~ 1 and the present day. The global SFR 
density is consistent with an evolution of the form (1 + z) 2 over < z < 1, with no 
evidence for a dependence on stellar mass. The sSFR of these star-forming galaxies at 
z ~ 0.75 does not depend upon the density of their local environment. Considering just 
high-density environments, the low-mass end of the sSFR-M* relation in our data is 
steeper than that in Stripe 82 at z = 0, and shallower than that measured by ROLES 
at z = 1. Evolution of low-mass galaxies in dense environments appears to be more 
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^ rapid than in the general field. 
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1 INTRODUCTION 

In recent years, evidence of the bimodal nature of the galaxy 
population has been obtained with increasing precision (e.g. 
Strateva et aLl|2001| |Baldry et aT||2004| |Bell et al.|[2007| ). 
Locally, the galaxy population divides quite cleanly into 
those which are actively star-forming and those in which 
star-formation has been terminated, or "quenched" . The rel- 
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ative mix of these two populations appears to be strongly 
dependent on both environment and stellar mass (M*, e.g. 
|Baldry et al.|2006| [Peng et al.|2010| ). In particular, the high 
mass end of the galaxy stellar mass function (GSMF) is 
dominated by passively evolving galaxies, while the actively 
star-forming popul ation dominates at st ellar masses below 
M - 10 10 M o (e.g. IPozzetti et al.||2010k at z<l. Observa- 
tions suggest that star lormation is truncated first in the 
most massive galaxies (e.g. Cowie et al.||1996 Bundy et al 



2006 



Pozzetti et al-poTO ); however, the stellar mass func 
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tion of actively star-forming galaxies itself evolves very little 
(e.g. |Pozzetti et aL|2010| |Gilbank et aT1|201Qa|b| ). 

Similarly, in the local universe the relative fraction 
of star-forming galaxies is strongly dependent on environ- 
ment, with the densest environments dominated by passive 
or "quenched" galaxies, and star- forming galaxies preferen- 
tially residing in lower density, "field" environments. But the 
properties of star-forming galaxies themselves have at most 
a weak dependence on environment (Balogh et al. 2004 Wolf 
|et al.|[2QQ9j |Vulcani et al.|[20To| . Recently, several authors 
have claimed evidence for evolution in this environment de- 



pendence (e.g 


. Gerke et al.||2007| |Cucciati et al.||2010 


|Bol- 


zonella et al. 


|2010| |McGee et al.||2011| |Patel et al.| 


201l| 


George et al. 


201lb, with several observations possibly in- 



dicating enhanced star formation rates (SFRs) in dense re- 
gions at z~l, under some circumstances (|Elbaz et al.[2007 



Cooper et al.||2008| |Ideue et al.||2009| |Sobral et afl|20TT 
Griitzbauch et al.||2011| ). All of these effects are relatively 



subtle, so comparisons between works are complicated by 
different definitions of environment (e.g. Cooper et al.| 2010 



Muldrew et al 



2012| ), sample selection and choice of esti- 
mator (e.g., average SFR, average specific star formation 
rate sSFR — SFR/M*, SFR of star- forming population), 
and star-formation indicators 



(e.g. 



Gilbank et al. 



2010a 



Patel et al. 2011). Indeed, previous, apparently contradic- 
tory, results may be reconciled when uniform definitions are 
ado pted (|Cooper et al.||2010| |Sobral et al~1|2011 ) . 

Peng et al.| ( 2010| ) recently presented an illuminat- 



ing, phenomenological description encapsulating the envi- 
ronmental and stellar mass dependence of galaxy activity 
and suggested that these effects appear to be entirely sep- 
arable. In their model, the efficiency of environment-driven 
transformation is independent of stellar mass and redshift, 
and the shape of the stellar mass function (SMF) for star- 
forming galaxies is universal and time-independent. How- 
ever, their model says nothing about the rate at which galax- 
ies transform from the star-forming to passive sequence; if 
this rate is slow enough, it will be observable as a popula- 
tion of primarily low-mass galaxies with lower-than-average 
SFR. 

A direct measurement of this timescale, which would 
provide important insight into the mechanisms driving this 
evolution, can be obtained by detecting a population of 
galaxies currently under the influence of "environment- 
quenching". The most likely place to find such a signature 
is amongst low-mass galaxies (for which mass-quenching is 
ineffective), at redshifts z > 0.5, when gas fractions and in- 
fall rates are high. Most spectroscopic surveys at these red- 



shifts are limited to fairly massive galaxies (e.g.|Noeske et al. 


2007 


Cooper et al.|2010 


|Bolzonella et al.|2010| |Patel et al. 


2011 


Muzzin et al.|2012 


). The Redshift One LDSS-3 Emis- 



sion Line Survey (hereafter ROLES) was designed to extend 
this work to lower stellar masses at z ~ 1, by searching 
for emission lines in K— selected samples, from fields with 
very deep imaging ( |Davies et al. 2009| . This was a spec- 
troscopic survey, conducted using the LDSS-3 instrument 
on the Magellan (Clay) telescope in Chile. With a custom 
made KG 750 filter, redshifts and [Oil] emission line fluxes 
were obtained for galaxies at 0.889 < z < 1.149 in the mass 
range 8.5 < log (M*/M ) < 9.5. 

ROLES demonstrated that the sSFR-mass relation 
evolves steadily with redshift, in a nearly mass-independent 



way, so the SFR density (SFRD) evolution is characterised 
primarily by an evolution in normalization only (Gilbank 
et al.|2010b ). However, there is a hint that the low-mass end 
of the sSFR-mass relation becomes steeper at z ~ 1 (| Gilbank 



et al.|20fT| ) , suggesting that the lowest-mass galaxies formed 



their stars later, and on longer timescales. Surprisingly, de- 



spite the small fields covered, |Li et al. (2011) found a clear 



environmental dependence amongst the star-forming pop- 
ulation. Star-forming galaxies in only moderately (factor 
~ 15) overdense regions at z = 1 appear to have higher 
SFR, a result that is opposite to the (weak) trend seen lo- 
cally. This is qualitatively consistent with results from some 
other surve ys (e.g. |Elbaz et al.|2007||Ideue et al.|2009||Sobral 
eTaLl[20TT| 



Here, we adopt the ROLES methodology (Gilb ank et al.| 



2010b) to explore star formation and its environmental de- 
pendence amongst low-mass galaxies over the redshift range 
0.62 < z < 0.885. This provides an intermediate link be- 
tween ROLES at z = 1 and the local Universe, using con- 
sistent galaxy selection and SFR measurement methods. 
Moreover, the redshift range and fields were chosen to in- 
clude highly overdense regions, including the well- studied 
MS1054-03 galaxy cluster (e.g. 



van Dokkum et al. 2000 



Forster Schreiber et al. 1120061). Thus the data span a wider 



range in environment compared with the ROLES data. 

This paper is presented as follows, ^describes the sur- 
vey and image reduction methodology, while details of the 
emission line detection procedure are presented in ^3] The 
basic measurements, corrections, and limiting values are pre- 
sented in Q Our results are shown in ^5] and we compare 
our results on the environmental independence of sSFR with 
published results at z = and z = 1 in ^6] Finally, our 
conclusions are summarized in ^7| AB magnitudes are used 
throughout unless otherwise stated and we use a ACDM 
cosmology of Ho = 70 km s _1 Mpc -1 , VLm — 0.3, and 
Q A = 0.7. Finally, note that all ROLES SFRs have been 
corrected using the empirical stellar mass dependent rela- 



tionship determined in |Gilbank et al. ( |2010a[ ), and described 
ing] 

2 DATA ACQUISITION & REDUCTION 

The design and implementation of the present survey is sim- 



ilar to our previous work at z = 1 (Gilbank et al. 2010b 
hereafter referred to as ROLES). In this section we review 
the target selection criteria, observation strategy, and image 
reduction steps. 



2.1 Target Selection 

Targets were selected based upon their i^-band magnitudes, 
22.5 < K < 24, and their photometric redshifts as provided 
by |Forster Schreiber et aT] ( |2006[ FIRES) and Mobasher 
& Dahlen ( |2009[ CDFS)~During the initial survey mask 
design phase, photometric redshifts were used to prioritize 
those targets which were expected to lie within our redshift 
range of 0.62 < z < 0.885, considering the probability dis- 
tribution of the photometric redshift. Galaxies with large 
photometric redshift uncertainties, or which were expected 
to lie outside our target redshift range, were also included in 
the mask design, with lower priority. As with ROLES, the 
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Figure 1. ROLES KG650 pointings in CDFS and FIRES. The ROLES pointings in CDFS are centered at 
(RA, Dec) = (03^32 m 27.6 s ,-27 d 45 m 00 s ) and (03^32 m 28.8 s , -27 d 52 m 12 s ) while the pointing in FIRES is centered at 
(10^56 m 58.26 s ,-03 d 37 m 0.53 s ). Each CDFS pointing is limited by the 8.2 arcminute diameter field of view (FOV) of the LDSS- 
3 spectrograph, shown as the thick black circle. The FIRES area is a 5.5 x 5.3 arcminute region which falls completely within the 
LDSS-3 FOV. All of the galaxies targeted are shown as open circles while those galaxies with observed emission lines (not necessarily 
[Oil]) are overlaid with filled circles. 



high sampling density in these fields means that the details 
of the prioritisation play a limited role, and the photometric 
redshift preselection does not constitute a strong prior. 

Our survey consists of two pointings in the Great Ob- 
servatories Origins Deep Survey (GOODS) region of the 
Chandra Deep Field South (CDFS, e.g. |Wuyts eral][2008 ) 
and one pointing in the MS 1054-03 cluster region of the 
Faint Infra- Red Extragalactic Survey (FIRES, e.g. |Forster 
Schreiber et al.|2006| jCrawford et al.|2011| ) field, for a total 
of 18 masks. The three pointings are shown in Figure^ with 
observational targets and emission line detections (described 
below) indicated. 



100 




6000 6200 6400 6600 6800 7000 
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2.2 Observations 

All spectroscopic observations were obtained using the 6.5 
meter Magellan (Clay) telescope. Multi-object spectroscopy 
for our 1946 targets was provided by the Low Dispersion Sur- 
vey Spectrograph 3 (LDSS-3). The spectra were dispersed by 
the medium red grism (300 lines/mm) which has a dispersion 
of approximately 2. 65 A/pixel at 6500A and a relatively uni- 
form throughput across the KG650 wavelength range. Com- 
bined with the plate scale of 0. 189" /pixel and survey mask 
slit width of 0.8", the resolution is 11.2A FWHM. 

The spectral wavelength range was limited to approxi- 
mately 650 =b 50 nm by a filter, herein referred to as KG650. 
The transmission curve for the KG650 filter is shown in Fig- 
ure^ From this transmission we define our sensitivity range 
as 6040A < X obs < 7025A. 



Figure 2. The transmission curve for the custom designed KG650 
filter. The vertical dashed lines at A = 6040, 7025 A indicate the 
spectral range considered by our survey, and correspond to the 
full width at half maximum (FWHM) locations of the transmis- 
sion curve. Since the survey specifically targets [Oil] emission at 
3727A, the survey redshift range is also restricted by the filter 
FWHM and is taken to be 0.62 < z < 0.885. 



The design of the survey masks was driven by the 
Nod fe Shuffle (NfeS, |Glazebrook fe Bland-Hawthorn|2001 



Gilbank et al.pOlOb ) observing strategy. The principle ad- 
vantage of this technique is that it allows for accurate sky 
subtraction at red wavelengths, where the sky brightness 
is dominated by rapidly varying emission lines. Target slits 
were 0.8" wide by 3.0" long, which allowed for nearly 200 ob- 
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Table 1. Spectroscopic masks with their corresponding target 
counts, total exposure times, and typical guide star seeing con- 
ditions. The two pointings in CDFS are labeled as CDFS.l and 
CDFS.2. The total number of unique survey targets is 1946. The 
number of targets, iVtargets, refers to the number of usable main 
survey spectra, excluding filler objects but including duplicates. 



Mask ID 


Field 


-^targets 


Exposure Time 
(hrs) 


Seeing 
(arcsec) 


mask24 


FIRES 


115 


2 


0.5 


mask25 


FIRES 


112 


2 


0.6 


mask26 


FIRES 


100 


2 


0.63 


mask27 


FIRES 


95 


2 


0.53 


mask28 


FIRES 


99 


2 


0.6 


mask29 


FIRES 


105 


2 


0.72 


mask30 


CDFS.l 


191 


2 


1.0 


mask31 


CDFS.l 


180 


2 


0.79 


mask32 


CDFS.l 


178 


2 


0.8 


mask33 


CDFS.l 


175 


2 


1.0 


mask34 


CDFS.l 


178 


2 


1.0 


mask35 


CDFS.l 


171 


2 


0.66 


mask36 


CDFS.2 


167 


2 


0.86 


mask37 


CDFS.2 


167 


2 


0.92 


mask38 


CDFS.2 


164 


2 


0.88 


mask39 


CDFS.2 


160 


2 


1.19 


mask40 


CDFS.2 


156 


2.5 


0.82 


mask41 


CDFS.2 


154 


2.5 


0.74 



jects to be targeted per mask given the LDSS-3 FOV of 8.2' 
diameter, and the requirement that dispersed object spectra 
must not overlap. We used a N&S cycle time of 60 seconds 
and a slit nod distance of 1.2". 

The exposure time for each mask is shown in Table [I] 
From the 18 masks combined, a total of 3292 objects were 
targeted. This includes objects which were targeted in multi- 
ple masks, and filler objects that do not satisfy our primary 
selection criteria. The total number of primary targets was 
2770 (including duplicates), of which 1946 are unique. These 
are split between 640 in FIRES and 1306 in CDFS. There 
were cases where slits were cut too close to mask edges, or 
to each other, and thus yielded unusable spectra. The num- 
ber of bad slits per mask typically ranged between 2 and 7. 
CDFS masks 40 and 41 had an unusually large number of 
bad slits, with 13 and 14 respectively. For the FIRES field, 
three exposures per mask were acquired while for CDFS 
we adopted a shorter exposure time per mask and acquired 
four frames per mask. The total exposure time per mask was 
identical for FIRES and CDFS except for CDFS masks 40 
and 41 where we acquired an extra 30 minute exposure for 
each due to time remaining in the observing programme. 



2.3 Image Reduction 

The FITS image files created by LDSS-3 were processed 
through an image reduction pipeline similar to that de- 
scribed in the Carnegie Observatories COSMOS (Carnegie 
Observatories System for MultiObject Spectroscopy) Cook- 
boolQ with custom- written IDL routines to supplement the 
existing software when required. 



See http:/ /obs. carnegiescience.edu/Code/cosmos/Cookbook. html 



2.3.1 Initial Frame Combining 

LDSS-3 was read out in two amplifier mode meaning that 
each mask exposure consisted of two raw FITS images. 
These images were combined using the COSMOS "stitch" 
routine with gain parameters set according to each specific 
amplifier and dewar parameter set to LDSS3-2. The "stitch" 
routine removes bias and corrects for differences between 
amplifier gains in LDSS-3 so no further bias removal was 
necessary after this stage. 



2.3.2 Bad Pixel Mask from Charge Traps & Cosmic Rays 

A bad pixel mask (BPM) was created using Nod & Shuffle 
dark frames (N&S darks) acquired using the same N&S pa- 
rameters as for the science frames. The resulting image is 
one which is mostly dark (predominantly read noise counts) 
with streaks of bright pixels indicating bad pixels created 
by charge being trapped by individual pixels in the CCD 
array. The streaks correspond to the charge shuffle distance 
and direction defined by the N&S strategy. A BPM was 
made from a N&S dark by dividing the original N&S dark 
frame by a 1x3 boxcar smoothed version (smoothing done 
in the direction perpendicular to the shuffle direction) of the 
same frame. Bad pixels appeared brighter in the ratio of the 
frames and were recorded as being bad in the BPM. Since 
several N&S darks had been acquired, they were each pro- 
cessed in the same manner and finally median combined into 
one single BPM. 

As many masks were observed with only three expo- 
sures, we used the IRAF COSMICRAYS task to identify the 
locations of cosmic rays (CR) in each frame based upon user- 
defined threshold levels and cosmic ray shapes. The pixel 
locations were recorded as a unique BPM for each exposure, 
and this resulted in superior quality image stacks during 
median combination. 



2.3.3 Wavelength Calibration 

The COSMOS apertures routine was used to make pre- 
dictions of initial positions of the slit centers in each mask. 
These positions were compared to the actual slit center po- 
sitions as imaged through the optical path and corrected 
(to less than 1 pixel difference) using the ALIGN-mask rou- 
tine. The positions of known arc lines were predicted for the 
arc calibration frames using the ALIGN-mask, map-spectra, 
and spectral-map routines. The COSMOS adjust-map 
routine was adequate for providing an initial wavelength 
calibration solution for most slits in a given mask. However 
there remained several cases where analysis of sky emission 
lines revealed inaccurate calibration. For this reason an IDL 
routine was used to determine a third order wavelength cali- 
bration solution to all of the slits in each mask, based on the 
position of these emission lines. Final emission line position 
residuals were typically < 0.7A. 



2.4 Creation of Stacked Frames 

For most masks the individual exposures to be stacked were 
acquired on different dates. As masks were interchanged in 
the optical path frequently and the telescope was at differ- 
ent orientations while tracking the target field at different 
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times of the year, differences in mask flexure, rotation, and 
shifts were introduced between one exposure and another. 
The transformations between science frames and an arbi- 
trary reference frame were determined based on the com- 
mon sky emission line centroid positions in each frame. The 
transformations commonly required a small rotation, shifts 
in the Xccd and Yccd directions, and on occasion a mul- 
tiplicative scaling. The IRAF task GEOMAP computed these 
transformations while the geotran task was used to apply 
them to each non-reference frame to be stacked. 

The applied geomap/geotran transformations ac- 
counted for differences in slit positions from one exposure 
frame to another. However, there were also cases where the 
target galaxy within a slit varied slightly in position be- 
tween the frames to be stacked. To rectify this, another IDL 
program was written to determine and apply any further 
required shifts in the spectral and spatial directions based 
upon a list of bright emission lines identified by eye and 
found in each frame to be stacked. Any such shifts were 
typically ~ 1 pixel in the Xccd and/or Yccd directions. 

The steps required to create the stacked signal frame 
from the N&S observations are described below. Individual 
exposure frames are labeled as A and B, and the recipe can 
be extended to an arbitrary number of frames. 

(i) Shift frame A by 16 pixels in the spatial ("y") direc- 
tion, to get a new frame A±; 

(ii) Perform the subtraction A - A\ to get a new frame 
A2. This is the sky subtracted frame; 

(iii) Shift frame A 2 by 6 pixels in the spatial ( "y" ) direc- 
tion, to get a new frame A3; 

(iv) Perform the subtraction A2 - A3 to get a new frame 
A4. This frame is the "positive" image frame; 

(v) Repeat steps (i) through (iv) for all individual expo- 
sure frames; 

(vi) Determine if individual frame flux scaling is neces- 
sary for each frame based upon the flux level ratios of sev- 
eral manually identified emission lines common between the 
brightest frame and the frame being scaled; 

(vii) Apply further (small) frame shifting if necessary 
based upon the centroided positions of the identified emis- 
sion lines used in step 6; 

(viii) Median add the "positive" image frames, A4 + B4 

a 

The IRAF task IMCOMBINE was used to median stack the 
individual "positive" image frames. The bad pixel mask de- 
scribed in §2.3.2} which includes identified cosmic rays, was 
used to ignore pixels during the combination. 

A corresponding noise frame was created in a manner 
similar to the stacked signal frame, as described by Gilbank 
et al. (2010b|. In the equations below the subscripts ij refer 



to the ij pixel of the frame. 

(i) Apply the same frame flux scaling, determined in step 

(vi) of the stacked science frame creation recipe, to each sky 
added frame; 

(ii) Apply the same (small) frame shifts, determined from 
the locations of common bright emission lines used in step 

(vii) of the stacked science frame creation recipe, to each sky 
added frame; 

(iii) Stack (median add) the sky added image frames to 
get a new frame, | < sky > |; 



(iv) Add in the LDSS-3 read noise, R. The read noise 
must be added in twice since the median combined frame 
consists of a shifted frame added to a non-shifted frame, 
each containing read noise. The read noise adjusted frame 
is calculated as follows: 

<sky > ^~) 2 + 2(i? 2 ) (1) 

(v) Scale frame Ni n div,ij by the number of individual sci- 
ence frames used in the median combination, n/ 

N in 



frames •> 



N c , 



(2) 



(vi) Shift the frame N CO m,ij by 6 pixels in the spatial ( "y" ) 
direction, to get a new frame N' corn ^j] 

(vii) Perform the quadrature addition of these last two 
frames to get the final noise frame: 



) 2 + (N' ..) 2 

/ 1 V com.ti J 



(3) 



3 EMISSION LINE DETECTION 

3.1 Creation of Signal-to-Noise Frame 

To identify faint emission lines, a normalized 2-D convolu- 
tion kernel, k ern , was created which had the same Gaussian 
shape as a typical bright emission line (FWHM= [5.5,3.5] 
pixels), and was convolved with the signal (S) and noise (N) 
frames to give flux- conserved, convolved signal and noise 
frames, according to: 



Sconv ,ij — Sij 



and 



A^" con 



Nl ® fcl„ 



(4) 



(5) 



The next step was to estimate the continuum found in 
the original signal frame. Similar to the convolution of the 
signal and noise frames, a convolution was again performed 
on the raw signal frame, using a 2-D normalized averaging 
kernel, k con t' 



Sij ® kcont • 



(6) 



The shape of the kernel consisted of a zero central region 
(20 pixels spectral by 3 pixels spatial) and two sidebands 
(also each 20 pixels spectral and 3 pixels spatial). The side- 
bands had the same Gaussian FWHM of 3.5 pixels in the 
spatial direction as the emission line kernel for their entire 
spectral length of 20 pixels. Convolving the kernel with the 
raw science frame provided an estimation of the continuum 
for the pixel located at the center of the kernel. The zero 
region was included so that the continuum estimate was not 
biased by the presence of an emission line. This provides a 
continuum estimate that is effectively an average of the flux 
in the spectral and spatial directions, in the "wings" of the 
pixel for which the continuum was being determined. 

The noise due to the continuum, N CO nt,ij, was calcu- 
lated by convolving the emission line kernel, k em , with the 
estimation of the continuum frame, C, as follows: 



cont,ij 



\/ Cij k^r, 



(7) 
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and the total noise frame, N to taU was calculated by adding 
in quadrature the convolved raw noise frame with the con- 
volved continuum noise estimate 



Ntota 



/V 2 4- TV 2 

%j ' cont,tj 



(8) 



In most cases the noise is dominated by sky line emission, 
but the continuum noise is not entirely negligible. 

This procedure accurately accounts for statistical noise 
in our spectra, but may not account for low-level systemat- 
ics resulting from weak charge traps, sky emission residuals 
(minimized but not entirely eliminated with N&S cycles of 
60s), or overlapping spectra. We analyzed the rms fluctua- 
tions in the final science frames and compared this with the 
associated noise estimate, for each mask. Specifically, thirty 
equally spaced "test" locations were chosen along the center 
line (distributed in the spectral direction) of each slit in a 
given mask, for both the stacked science and stacked noise 
frames. The mean pixel value for each slit test location in 
the science frame was determined by taking the mean of the 
pixel values within two 60 pixel sidebands, located to either 
side of the test location. The fluctuation of the test loca- 
tion pixel value from the mean was then simply the actual 
pixel value subtract the mean value. For every slit in the 
mask, the science fluctuation (as) and noise value (/i n ) for 
thirty test locations were recorded. Histograms of the science 
fluctuations and corresponding noise values were then fitted 
with Gaussians. Finally, the ratio of the best-fit Gaussian 
standard deviation of the science frame fluctuations and the 
Gaussian mean of the noise values gave the "noise correction 
factor" (NCF): 



NCF : 



(9) 



A typical noise correction factor was ^1.2, indicating that 
residual systematics amount to an additional 20% on top of 
the statistical noise. 

The final S/N frame was calculated as 



' S_ 
N 



■ Co. 



N, 



total, ij 



■NCF 



(10) 



3.2 Emission Line Finding 

The central five rows of each spectrum was extracted from 
the 2D frame, to minimize effects near slit edges that affect 
line detection. For every pixel above a S/N threshold of 3, 
an "n-connected neighbour" search was performed to locate 
connected neighbouring pixels that also exceed this thresh- 
old. A candidate detection then consists of two or more con- 
nected pixels; if multiple detections were separated by five 
pixels or less, they are combined into a single detection. De- 
tections found within three pixels of the spectral ends of the 
extracted spectrum, and those that were due to overlaps of 
th order spectra, were excluded. The resulting list was visu- 
ally inspected, and obvious false detections (due in general 
to overlapping spectra or missed cosmic rays) were manually 
removed. 



3.2.1 Catalogue Purity 

The 3a (S/N ^ 3) catalogue was internally tested in 
two ways to determine the purity, following Gilb ank et al.| 
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Figure 3. Main Panel: The cumulative fraction of the number 
of recovered (solid green curve) vs. spurious (dashed red curve) 
emission lines for galaxies which were targeted in multiple masks, 
as a function of significance. Inset: Histogram of the number 
of recovered and spurious emission lines. This demonstrates that 
95% of spurious detections occur below 5cr (highlighted as the 
solid vertical blue line). 
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Figure 4. Main Panel: The cumulative fraction of the number 
of recovered (solid green curve) vs. spurious (dashed red curve) 
emission lines for galaxies which contain at least two emission 
lines and lead to self-consistent redshifts, as a function of signifi- 
cance. Inset: Histogram of the number of recovered and spurious 
emission lines providing self-consistent redshifts for the galaxies 
containing them. It is clear from the cumulative fractions that 
~ 97% of spurious detections occur below 5<j (highlighted as the 
solid vertical blue line). 



(2010b). We first consider the reproducibility of emission 



lines for the 412 galaxies that were targeted on more than 
one mask. For these galaxies, detection lists were compared 
and emission lines were considered to match if their wave- 
lengths were within ±6.5A (2.5 pixels) of each other. If a 
detection was found in all of the masks the galaxy was tar- 
geted in, then it was considered fully recovered; otherwise 
it was considered spurious. This is therefore a conservative 
estimate of the purity. The results of this test are shown in 
Figure [3] where it is clear that 95% of spurious detections 
occurred in detections below 5cr. 

An independent test of purity is to consider spec- 
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tra for which more than one candidate is detected. The 
wavelengths of these candidates were compared with ex- 
pected sets of lines, which are likely to appear in only a 
small number of combinations: the Hp — [OIII] complex 

(Hp, [0///] 4 959, [Om] 5 007), and the [NeIII} 38 69, [OII} 3 727 

pair. Candidates were considered real detections if their line 
ratios matched one of these combinations. In the case that 
the lines did not correspond to an expected set, the line 
with the highest significance was considered real, while the 
next-highest significance line was considered spurious; lower 
significance lines were omitted for the purpose of this test. 
Figure [4] shows the S/N distribution of these real and spuri- 
ous lines; ~ 97% of spurious lines have significance less than 
5a. 

From these tests, we conclude that > 95% of false detec- 
tions occur below a significance threshold of S/N ^ 5cr; thus 
we only consider detections above this limit in our analysis. 

3.3 Emission Line Flux Determination 

Flux calibration was based on the spectrophotometric stan- 
dard star HD 49798 ( |Bohlin fe Lindler 11992) . Emission line 
fluxes and their errors were measured from the stacked raw 
science frames with the continuum estimation removed, and 
the stacked raw noise frames. The bad pixel masks were in- 
corporated to eliminate bad pixels and cosmic rays. For each 
detection, the centroid position of the flux was found within 
a 15 x 17 pixel box, initially centered on the location of 
the highest significance pixel in the emission line. The total 
emission line flux was taken to be the sum of flux within a 7 
x 5 pixel region about this centroid. The line flux error was 
calculated for the same pixels, based on the noise spectrum. 

To account for varying photometric conditions, we com- 
pare the flux in the continuum measured from the spectra 
on each mask with photometric data from public catalogues. 
For the CDFS field we use the R-band magnitudes from 
FIREWORKS flWuyts et al.|2008) , which covers almost ex- 
actly the same wavelength range as our spectroscopy. For 
the FIRES field, the available photometry does not include 
R-magnitudes, so we interpolated between the HST WFPC2 
F606W and F814W filters. For each mask we calculate the 
average offset between the flux in our spectra and the con- 
tinuum flux measured from the imaging. We use this to iden- 
tify the most photometric mask in each field, and the offsets 
from this mask for all of the others. We then correct the 
flux calibration for the non-photometric masks, to match 
this reference frame. The correction is typically ^0.5 mag, 
with a maximum of 1 mag. 

We take advantage of galaxies within our 5cr linelist 
that were imaged on multiple masks, to further check the 
consistency of the flux calibration and our uncertainty es- 
timates. The flux differences for separate observations were 
determined and plotted as a function of line flux, as shown 
in Figure [5] As expected, the matching line flux differences 
are scattered about zero. The significance of each difference 
is obtained by dividing by the flux uncertainties added in 
quadrature. The significance distribution has a standard de- 
viation of a ~ 1.28, but a Kolmogorov-Smirnov test cannot 
distinguish between this distribution and a normal distri- 
bution with a = 1. Thus we conclude that uncertainties in 
relative flux calibration from mask-to-mask are negligible. 

We looked for a correlation between the flux difference 
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Figure 5. The difference in ^ 5cr emission line fluxes for match- 
ing detections found in multiply-imaged galaxies as a function of 
line flux. The zero difference line is shown as the horizontal green 
dashed line. 



compared with the photometry and galaxy size. While there 
does appear to be a correlation in the expected sense, that 
larger galaxies are missing more flux in our spectroscopy, 
there is a lot of scatter from object to object. We have elected 
therefore not to apply an aperture correction, but note that 
the fluxes for our largest galaxies are likely underestimated. 
Only three of our spectra overlap with spectra obtained by 
Vanzella et al. (2008); while this comparison shows our flux 



calibration is consistent with theirs, there are not enough 
objects in common to state this with a high degree of confi- 
dence. Thus we expect our line fluxes are dominated by this 
systematic uncertainty in flux calibration, which is likely at 
least a factor of ^ 2 with a dependence on galaxy size. 

The wavelength- dependent flux limit was determined 
for each survey mask from the associated noise propagated 
through the analysis pipeline. For each mask, the average 
noise spectrum a\ was determined from all the dispersed 
spectra in the mask, and secure detections were then defined 
as those brighter than 5<ta . Figure [5] shows the flux of all 
detected lines, and the average ba\ noise level for all masks. 
Unless otherwise stated, analysis in this paper excludes lines 
that fall below the average flux limit shown here, to enforce 
a uniform limit. 



3.4 Line Identification 

For galaxies in the ^ 5cr catalogue where more than 
one emission line was detected with an appropriate wave- 
length separation, identification was straightforward. How- 
ever, there are only 41 such candidates, 15 of which are 
identified as [Oil]. The remainder of the catalogue consists 
of single emission lines, for which we rely on the photometric 
redshift probability distribution functions (PDFs) to iden- 
tify the line, as described in Gilbank et al.j ( |2010b ). The 



relative likelihood of a line being [Oil] was assigned to each 
detection in our 5cr catalogue by determining the ratio of 
the probability of the emission line being [Oil] to the to- 
tal probability of being either [Mgll], [CIV] or one of the 
Hf3 — [OIII] complex. The probabilities were calculated by 
integrating the photometric redshift PDFs over the redshift 
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Figure 6. Fluxes for emission lines detected with 5cr significance 
are shown as a function of wavelength. Detections are divided 
into the CDFS ROLES (solid/open red circles), FORS2 (blue 
diamonds), and FIRES (solid/open green squares) fields. The av- 
erage 5Nij noise flux limit of all masks is overplotted as a solid 
line. 
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Figure 7. The published redshift, where available, is shown for 
galaxies with candidate emission lines in our catalogue as a func- 
tion of the line wavelength. Angled lines show the relation ex- 
pected for seven different emission lines as indicated; in particu- 
lar the black dashed line indicates our target line of [Oil]. Filled 
circles correspond to detections with C[OII] > 0.9. 



ranges corresponding to the lines of interest in our spectral 
window. 

For a few lines, publicly available spectroscopy covering 
a larger wavelength range are available; nine from|Crawford| 
|et al 



J20TT| MS1054-03) and 73 from |Wuyts et al.| |2Q08 
FIREWORKS) have secure redshifts. Figure [7] shows these 
secure redshifts as a function of the wavelength of the de- 
tected emission line in our sample. Straight lines indicate 
the redshift- wavelength relation for different emission lines, 
including [Oil] (black dashed line). The detections denoted 
by solid symbols were those which had an integrated pho- 
tometric redshift probability of being [Oil]. Most of the 
emission lines which had a high probability of being [Oil] 
(C[OII] > 0.9) are confirmed as such from the public red- 
shift. Of the 82 matched detections shown in the figure only 



three appear inconsistent with their published redshift, sug- 
gesting that the purity of our sample is over 95 per cent. 

3.5 Final 5cr Catalogue 

The final catalogue contains all detections with S/N ^ 5, 
and a redshift determined following the line identification 
procedure described above. The detections were also given 
a line quality flag as follows: 

0: Line is not likely [Oil] based on the photometric red- 
shift PDF, and there is no existing public redshift 

1: Photometric redshift is consistent with the detection 
being [Oil] 

2: Photometric redshift is consistent with the detection 
being [Oil], confirmed by detection of [Nelll]. 

3: Photometric redshift is consistent with the detection 
being [Oil], confirmed by a published redshift. 



4 ANALYSIS 

Spectroscopic redshifts determined for the galaxies in our 
5cr line list, as well as photometry in the CDFS field from 



the FIREWORKS catalogue flWuyts et aL|2008| C7 38 , #435, 
B, V, y 60 6, R, ir?5, I, Z8 50, J, H, K a , [3.6], [4 5], [8.0]), 
and in the FIRES field ( |Forster Schreiber et aL||2006| U, 
B, V, V606, ^814, J s , H, and K 8 ), served as inputs to the 



stellar PEGASE.2 ( |Fioc fc Rocca-Volme range 1997k popu- 
lation models, as described in Glazebrook et al. (2004). The 



models were fit to aperture magnitudes, and the final stellar 
masse^Jwere subsequently scaled according to the total K s 
magnitude. 

Emission line luminosities were first converted to "fidu- 



cial" star formation rates (SFRf), starting from the Kenni- 
cutt| dl998| ) relation 

1 q0AA Hi 



SFR f (M Q /yr) 



0.5 



7.9 x 10" 
1.82 



[L [0 n](ergs/s)} 
(11) 



where the factor of 1.82 accounted for the conversion from 



a Salpeter IMF flSalpeter| [l955) to the BG03 IMF flBaldry 
fe Glazebr ook 2003), and we assume Ah a = 1. In practice, 
Ah a and the other coefficients in Equation [TT] are likely to 
depend strongly on the galaxy stellar mass. We therefore use 



the empirical correction advocated by |Gilbank et al. (2010a 
Eq.8), based on analysis of the SDSS: 



SFRiM^yr- 1 ) 



SFR f 



{(a) -tanh[(X + 6)/c] 



(12) 



where X = log (M*/M ), a = -1.424, b = -9.827, c = 
0.572, and d = 1.700. We apply this correction to all SFR 
reported in this paper. 

This SFR estimate assumes that the [O II] emission 
arises from gas ionized by massive stars; any contribution 
from an active galactic nucleus (AGN) would reduce the 
SFR. We would not expect AGN to contribute significantly 
in such low-mass galaxies; in Gilbank et al. (2010b) we con- 



firmed from analysis of mid-infrared colours and Chandra 

2 PEGASE.2 only accounts for luminous stars in its determina- 
tion of stellar mass. Thus stellar remnants such as white dwarfs, 
neutron stars, and black holes were not included. 



© 2011 RAS, MNRAS 000,[T}[T6] 



The evolution of star formation in dwarf galaxies 9 




0.0 0.5 1.0 

Veoe " '775 (mags) 



1.5 



Figure 8. The colour-magnitude diagram of the CDFS and 
FORS2 targets is shown. Upper Panel: Our targets in the CDFS 
are shown as solid black circles while the public FORS2 data are 
represented by black asterisks. Bottom Panel: The overall dis- 
tribution of V606 — ^775 colour is shown. The red sequence cut 
is shown as the vertical green dashed line; we exclude galaxies 
redder than this limit from the remainder of our analysis. 



point sources that AGN were entirely negligible in that z = 1 
sample. We will therefore also neglect AGN contamination, 
except when we include higher mass galaxies from other sur- 
veys (§141). 
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Figure 9. Upper Panel: Summation of all photometric red- 
shift PDFs within the LDSS-3 FOV and with 22.5 < K < 24 
(black curve) compared with the summation of those photometric 
redshift PDFs corresponding to galaxies targeted in this survey 
(lower green curve), for the CDFS field. Lower Panel: Same as 
the upper panel, corresponding to the FIRES field. 



are in the blue cloud, a subset of the FORS2 galaxies are on 
the red sequence, with V<506 — ^775 > 1-25 mag. 

For each FORS2 emission line we calculated Vmax as 
before (with a typical value of 4.7 x 10 4 Mpc 3 ). The K- 
magnitude binned FORS2 completeness was determined in 
the same way as for our data, and the results are shown in 
Figure [l0| We extracted spectroscopic redshifts, line identi- 
fications and quality flags, and ID spectra for these galaxies. 
[Oil] emission line fluxes were measured from the ID spectra 
in the same way as for our own data. A constant 4cr noise 
flux limit of 6 x 10 ~ 18 erg s -1 cm -2 was adopted, approxi- 
mately the same as the average noise flux limit for the rest 
of our sample (see Figure |6j). 



4.1 Extension of the survey to higher masses 



We supplement our survey with brighter (more massive) 
galaxies, primarily from publicly available VLT/FORS2 
spectroscopy overlapping our CDFS sample area, from 
(Vanzella et al. 2008). Their sample was a colour and pho- 
tometric redshift selected catalogue with targets found be- 
tween the redshift ranges of 0.5 < z < 2 and 3.5 < z < 6.3. 
We selected only those targets which were found within the 
LDSS-3 field-of-view centered on our CDFS field pointings, 
and which fell within our redshift range of 0.62 < z < 0.885. 
Their observation masks used \" slits (compared to 0.8" for 
ROLES) and exposure times for each mask were typically 
^ 4 hours. Hereafter this higher mass sample is referred to 
as FORS2. 

For these massive galaxies, there is more concern that 
the emission lines could arise from AGN or LINER emis- 
sion, rather than star formation. As in ROLES, we there- 
fore exclude red-sequence galaxies from the sample, using a 
colour-magnitude diagram (CMD) consisting of bands which 
bracket the 4000 A break at z ~ 0.75. Figure ^ shows the 
characteristic bimodal colour distribution of galaxies in our 
CDFS and FORS2 samples; while all our low-mass targets 



4.2 Completeness 

As for ROLES, we characterize our spectroscopic complete- 
ness as follows. For each field, all photometric redshift proba- 
bility distribution functions corresponding to galaxies within 
our target fields were first binned by if -magnitude. Within 
each bin the PDFs were summed, giving the total redshift 
distribution for all galaxies in each bin, Pk(z). The summed 
redshift distribution in each bin was then integrated over the 
redshift of interest here, 0.62 < z < 0.885. The process was 
then repeated for only those galaxies that were successfully 
targeted (i.e., the slit was successfully extracted), and the 
ratio of this to the total distribution yields the targeting 
completeness. The redshift PDFs of targeted and candidate 
galaxies are shown in Figure [5] 

The resulting completeness is high, ~ 70 per cent, and 
independent of K magnitude. This is shown in Figure [To] 
with the CDFS field represented by the red dashed line and 
the FIRES field denoted by the green dashed line. The figure 
also includes spectra for brighter galaxies from the public 
domain, discussed in § |4.1| 
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Figure 10. Completeness histograms for our surveys in the CDFS 
and FIRES are shown as the red and green lines. We also include 
existing FORS2 data for brighter galaxies, shown as the blue his- 
togram. 



4.3 Redshift distribution and definition of dense 
environments 

The redshift distribution of our emission line galaxies is 
shown in Figure Two prominent peaks in this distri- 
bution at z rsj 0.68 and z ~ 0.73 correspond to the well- 
known "wall-like" structures in the CDFS (Gilli et al.|2003 



Le Fevre et al.||2004| |Vanzella et al.||2005| |Ravikumar et al.| 
2007) . We associate all CDFS galaxies with \z - 0.668 1 < 
0.016 and \z — 0.735 1 < 0.009 with these structures. As traced 
by our low-mass galaxy sample, this structure is spread fairly 
uniformly over the LDSS3 field of view, without an appar- 
ent density gradient or central concentration. The rest-frame 
velocity dispersion of emission line galaxies in each of these 
structures is 970km/s and 430km/s, respectively. The other 
important structure, at z ~ 0.83, is the MS1054-03 clus- 
ter, in the FIRES field flForster Schreiber et aT1|2006| ); all 
galaxies in the field and within Az = 0.02 of this redshift 
are associated with the cluster. The rest-frame velocity dis- 
persion of these emission line galaxies is 1300km/s, in good 
agreement with |van Dokkum et al.| ( |2000 ). Together, these 
three subsets of galaxies are referred to as "dense environ- 
ments" for the subsequent analysis. Combined, the subsam- 
ple comprises 112 galaxies, 23 of which are associated with 
the MS1054-03 cluster. The remaining galaxies are referred 
to as the "field"; it now represents an underdense sample 
relative to the average of our full sample. We will show in 
§ [5T2l that both the CDFS "walls" and the MS1054 cluster 
are comparably overdense, by a factor of at least 7 relative 
to the field, and probably more like a factor > 45. 




0.6 0.7 OS 0.9 

Redshift 

Figure 11. The redshift distribution of emission line galaxies in 
our sample. The distributions are normalized to an area of 100 
arcmin 2 , close to the area of the CDFS and ~ 3.6 times larger 
than the FIRES field. Large-scale structure, due to the known 
walls in the CDFS and the cluster MS1054-03 in the FIRES field, 
is identifiable as three narrow redshift peaks. 
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Figure 12. K-magnitude as a function of stellar mass for emis- 
sion line galaxies in our two survey fields (red for CDFS, and 
green for FIRES), supplemented with the more massive FORS2 
galaxies (blue triangles). The solid lines show the K-magnitude 
and adopted stellar mass limits, while the blue line is the best fit 
to the three populations combined. 



4.4 [Oil] Luminosity, Stellar mass and SFR Limits 

The stellar mass limit of our sample is determined from the 
scatter in the correlation between X-magnitude and stellar 
mass shown in Figure [l2| The horizontal line shows our lim- 
iting selection magnitude of K = 24. Based on the scatter 
in this relation, the sample is nearly (> 90 per cent) com- 
plete in stellar mass for log (M*/M©) > 8.85. We take this 
to be our 2a mass completeness limit; the sample extends 



to lower masses, but is systematically missing galaxies with 
high M/Lk ratios. 

The average 5cr [Oil] flux limit as a function of wave- 
length was shown in Figure [6] The sample is statistically 
complete for fluxes as low as ~ 5 x 10 _18 ergs s _1 cm -2 . Con- 
sidering our low redshift bound of z ~ 0.62, the correspond- 
ing [Oil] luminosity limit at which the sample should be 
statistically complete is determined to be log L[qii] ~ 39.9. 
However, most of the volume is limited to higher luminosi- 
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Figure 13. The SFR of all emission-line galaxies in our sample 
is shown as a function of their stellar mass. Our LDSS3 data are 
shown as the red points; blue points are public data from FORS2 
spectroscopy. The 2a stellar mass and 5cr SFR limits are shown 
as solid, black lines. 



ties, and log L[ ii] ~ 40.1 is a more representative limit 
for most of the data. The limiting SFR (including the mass- 
dependent empirical correction) can be determined from the 
[Oil] luminosity limit, using Equations 



11 



and 12 Figure 13 



is a plot of the empirically corrected SFRs versus stellar 
mass. The 2a mass and 5cr SFR limits are indicated with 
solid lines. We reach SFR > 0.1 Mq at the low stellar masses 
of interest here, corresponding to a mass doubling time (as- 
suming a recycling factor R = 0.5) of td °° 1.3 x 10 10 years. 
This is almost twice the Hubble time at z = 0.7, and we 
expect this depth is sufficient to capture most of the star 
formation at these masses (e.g. |Noeske et aT|2007| |Gilbank| 
etaLl[20TT| . 



4.5 Survey Volume and density estimates 

Our survey volume is denned by the survey area, the lim- 
iting magnitudes (22.5 < K < 24), and a flux limit on the 
emission line detection. As in ROLES, we determine the 
redshift limit for which the galaxy would fall outside our K 
magnitude limits, including the /c-correction term 



-2.58z + 6.67z 2 - 5.73z 3 - 0.42z 4 



(13) 



1 - 2.36z + 3.82z 2 - 3.53z 3 + 3.35z 4 ' 
The wavelength-dependent flux limit shown in Figure [6] 
was similarly used to determine the redshift limits over 
which each detected emission line would be observable. The 
volume, Vmax, from which a galaxy with a detected emission 
line could have been found was then calculated from the sur- 
vey area (105.62 and 29.15 square arcminutes for our CDFS 
and FIRES pointings, respectively) and the redshift space 
bounded by the K-magnitude and noise flux limits, and the 
wavelength limits of our spectra. The volume, V m ax i was 
determined for each galaxy by integrating the differential 
co- moving volume (see Hogg 1999} between the appropriate 
redshift limits. The maximum (0.62 < z < 0.885) Vmax is 
4.7 x 10 4 Mpc 3 in the CDFS fields, and 1.3 x 10 4 Mpc 3 
in FIRES. 

The three structures are defined by the redshift limits 
given in § |4.3| Interpreting the redshift limits as cosmologi- 
cal, removing these regions reduces the volume of our field 



sample by ~ 20 per cent. For the MS1054 cluster and CDFS 
walls themselves, which are decoupled from the Hubble flow, 
most of the galaxies are likely located within a volume that 
is much smaller than this cosmological volume. We will as- 
sume their line-of-sight extent is 10 Mpc (i.e. assuming a 
~ 5 Mpc virial radius, which is still probably conservatively 
large). This corresponds to a volume of ~ 200Mpc -3 for 
the cluster, and ~ 520Mpc~ 3 and ~ 610Mpc~ 3 for the two 
CDFS walls. 



4.6 ROLES and SDSS Stripe 82 data 

We will compare our results with similarly-selected data, 
from ROLES flGilbank et al.|2010b) at z ~ 1, and the local 



Universe from our Stripe82 analysis (Gilbank et al. 2011). 



Both of these samples are consistent with our present anal- 
ysis in the choice of IMF, the empirical calibration of [Oil] 
to SFR, and the removal of red, massive galaxie^] 

The ROLES SFR are also limited by [O II] flux, and 
thus the mass-dependent SFR limit has the same form as 
shown in Figure |13| However, the greater luminosity dis- 
tance and brighter sky at the wavelength of redshifted [O II] 
at z = 1 means, despite the longer exposure times, that the 
limiting SFR at z = 1 is about a factor ~ 3 greater than in 



the present study. In I Gilbank et al. (2011 ), we demonstrated 



that, locally, the SFRD has converged for SFR > 0.1M© /yr, 
for galaxies with log (M/M ) ~ 9, and SFR > lM /yr for 
galaxies with log (M) ~ 10. Globally, the average sSFR is 
known to evolve approximately as (1 + z) 2 ' 5 (e.g. |Prescott 
et al. 2009). Assuming the shape of the SFRD does not 



evolve strongly, we would expect ROLES (log (M/M ) > 9 
at z = 1) to be complete at ~ O.5M /yr, and the present 



study (log (M/M ) > 9 at 



0.75) to be complete at 



SFR rsj O.4M /yr. Thus, we expect our samples are deep 
enough to have recovered most of the star formation in the 
Universe, and to be fairly insensitive to the precise choice of 
limiting SFR. 



5 RESULTS 

5.1 Star formation rate density 

The SFRD was computed as follows: 



Psfr(M*) = 



Pou,i ■ SFRj 

ax, % 



(14) 



where the SFR was calculated according to Equations[TT]and 
|12[ and the sum is over all galaxies in a given bin of stellar 
mass. Pon,i is the probability that the line is [Oil], relative 
to the probability of it being any other plausible emission 
line; w% is a magnitude-dependent weight to account for in- 
completeness. Figure ^] shows the SFRD for our LDSS3 
sample. Recall from Figure [l0| that our completeness drops 
significantly for bright (K < 22.5) galaxies, corresponding 
to M > 3 x 1O 9 M . Moreover, this high mass end has no 



3 In fact, the Stripe 82 analysis is based on Ha measurements 
of the SFR which, by construction, are consistent with the mass- 
dependent [O II] -SFR conversion. This avoids the incompleteness 
in [O II] noted by Gilbank et al. (2011 ) for high mass galaxies in 
Stripe 82. 
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Figure 14. The SFRD of our sample, at z ~ 0.75, is shown 
as the red points with error bars at M < 3 X 1O 9 M0. The grey 
vertical dashed line highlights the mass limit of the present survey. 
The high-mass data are taken from the stacked radio analysis of 
galaxies at 0.6 < z < 0.8 from |Karim et al. ( |2011| ), shown as 
the open circles. Note that |Karim et al.| ( |2011 > find a factor ~ 2 
difference in normalization between their radio-based SFRD and 
the [OII]-based analysis of ROLES at z = 1. If this is treated as 
a systematic effect, the solid circles of |Karim et al.| ( |2011| on this 
Figure should be brought down by a factor two. This SFRD is 
compared with ROLES at z = 1 (blue diamonds), and the SDSS 
from our analysis of Stripe82 data (green triangles). The solid 
and dashed line represent Schechter functions fit to the Stripe82 
data, scaled by a factor (1 + z) 2 to z = 0.75 and z = 1. 



contribution from the FIRES field, and is therefore espe- 
cially sensitive to cosmic variance (not included in the er- 
ror bars) given the large scale structure in the CDFS field. 
Finally, no aperture corrections have been applied, which 
means SFRs are likely underestimated for the most massive 
galaxies. Thus we limit our data to M < 3 x 1O 9 M , and re- 
place the higher mass measurements with the 0.6 < z < 0.8 
data of |Karim et al.| ( |2011| ), who used stacked radio obser- 
vations in the COSMOS field to measure SFR over a wider 
field than we have here. They note that there is a factor ~ 2 
difference in normalization between [OII]-based SFR and 
their analysis. To eliminate this likely systematic offset, the 
Karim et ah] ( |2011| ) SFRD should be reduced (or our [Oil] 
data should be increased) by this factor. However, we show 
the data as measured, with no rescaling applied. 



We compare this with results from ROLES (Gilbank 
|et al.|[2010b| at z ~ 1, and the local Universe from our 
Stripe82 analysis (Gilbank et al. 2011| . Interestingly, the 
ROLES SFRD is consistent with our new data at z ~ 0.75, 
over most of the mass range, with some evidence for higher 



SFR at M > 1O 1U M if the SFRD from Karim et al 
are reduced by a factor 
with the z 



(2011) 



2 as described above. Compared 
0.1 data from Stripe 82, the SFRD at both 
z = 0.75 and z = 1 is higher by a factor ^ 3 at all 
masses. Given the statistical and systematic uncertainties, 
though, we cannot rule out a continuous evolution of the 
form SFRD oc (1 + z) 2 ' (e.g. |Gilbank et al.||2010b| ), repre- 
sented on Figure [l4| as the solid and dashed lines at z — 0.75 
and z = 1, respectively. 

Our findings are in good agreement with |Karim et ah] 



(2011 ) and |Gilbank et aL ( 2010b| ), that there is no evidence 
here for strong evolution in the shape of the SFRD. We note 



that both ROLES and the present work hint at an increased 
contribution from the lowest-mass galaxies for which we are 
complete (~ 1O 9 M ), relative to the local Universe. How- 
ever, any mass-dependent effect is fairly subtle, and it is un- 
likely that systematic uncertainties are understood to this 
level across the entire mass and redshift range. The data are 
not inconsistent with a simple, mass-independent evolution 
of SFRD oc (1 + z) 2 - over the redshift range < z < 1, as 



suggested by |Gilbank et~aT| ( |2010b ) 



5.2 Stellar mass function of star-forming galaxies 

The galaxy stellar mass function for our sample of star- 



forming galaxies (SF-GSMF) is shown in Figure 15 In or- 
der to study its evolution, we compare our z^0.75 data with 
the Stripe 82 SDSS data as detailed in § |4.6| The left panel 
only considers the present z^0.75 data above a conserva- 
tive SFR limit of 0.3 M yr -1 so as to compare directly 
with the depth of the slightly shallower z^l ROLES data. 
To compare fairly with the local data, the Stripe 82 sample 
is limited to SFR ^ 0.1 M yr" 1 . Since the SFRD falls 
by a factor of ~3 over this redshift interval (as discussed 
in § 5.1), this cut corresponds to a similarly evolving limit. 



For comparison, we also show the single Schechter function 
fit to the SDSS star-forming population at z = 0.1 from 



Peng et al. ( 2010 ) . Transforming to the |Baldry &; Glazebrook 
(2003J) IMF, the parameters of this fit are logM* = 10.92 
and <£* = 2.612 x 10~ 3 dex _1 Mpc" 3 . The high-mass end of 
the SF-GSMF function (M>3x 1O 9 M ) is not very well- 
determined in our present data, due to the limited survey 
area and spectroscopic sampling. At lower masses, where we 
expect our sample to be highly complete, the data are con- 
sistent with no evolution of the SF-GSMF from < z < 1, 



in agreement with the conclusions of Gilbank et al. (2010b) 



and Peng et al. (2010) 



0.75 

i 



In the right panel, the full depth of the present z 
dataset is considered, down to a SFR limit of O.1M yr 
The Stripe 82 data from the left panel is reproduced here, 
but we also show these data extended to a lower limiting 
SFR^O.O3M yr" 1 , allowing for a fa ctor ^3 in SFR evolu- 
tion. We also compare with data from |Pozzetti et aT| ( |2010| ), 
shown without error bars, for clarity. This sample defined 
star-forming galaxies as those with log sSFR/Gyr > — 1, 
approximately consistent with the limits we apply to our 
data, here. Again we convert the result to correspond to 
a BG03 ( |Baldry fe Glazebrook||2003| ) IMF. The results are 
consistent with little or no evolution in the SF-GSMF over 
this redshift range. However, it also demonstrates the sen- 
sitivity of any measured evolution to the limiting SFR of 
the samples. If we account for the global evolution of SFR 
we find the SF-GSMF remains constant down to the low- 
est masses for which we have statistically complete samples. 
Choosing a fixed, non-evolving limit would result in a large 
decrease with increasing cosmic time of the SF-GSMF at the 



low mass end. The difference between the analysis of [Peng 
et al. (2010), where SF galaxies are identified strictly by 



colour, and our Stripe82 analysis, shows that remaining sys- 
tematic uncertainties of this type are at least as important 
as any physical evolution. 

In Figure [16] we divide our sample into different en- 
vironments: the MS 1054 cluster in FIRES, the large-scale 
structure in CDFS, and the remaining population which we 
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Figure 15. Left panel: The stellar mass function of all star-forming galaxies with SFR^0.33Mq yr — 1 in our sample is shown as the 
points with error bars. This is compared with ROLES data at z ~ 1 (blue asterisks), to the same SFR limit, and with the SDSS Stripe 82 
data to a limit of a factor ~ 3 lower to account for evolution (see Gilbank et al. 2010b). The dashed line represents the Schechter function 
fit to the local, SDSS data byK ng et al.| ( |2010| . Right panel: Our data are now shown to a deeper SFR limit of SFR^O.IM© yr 1 ; this 
is compared to the Stripe82 data to two different depths, as described in the legend. We also compare with data at 0.55 < z < 0.75 and 
0.75 < z < 1.0 from |Pozzetti et ah] (2010), shown as black and red crosses respectively. The thick grey vertical dashed line in both panels 
highlights the mass limit of the present survey. 
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Figure 16. The stellar mass function for star-forming galaxies 
in our survey is shown, divided by environment as shown in the 
legend. Here we assume that the line-of-sight extent of both the 
cluster and the CDFS structures is 10 Mpc (comoving). The lo- 
cal Schechter function for star-forming galaxies from |Peng et al.| 
( |2010| ) has been renormalized to fit the low-density and struc- 
tured environments. We see no dependence of the shape on envi- 
ronment. Both the cluster and the CDFS "walls" are comparable 
in overdensity, a factor ~ 45 times denser than the field. Error 
bars and upper limits are la. The lower limit on the most mas- 
sive FIRES point is due to the lack of redshifts for K < 22. The 
vertical dotted line highlights the mass limit of the present survey. 



call the "field". The cluster in FIRES has no public spec- 
troscopy for K < 22, which me ans that field is incomplete 
for M > 10 9 - 5 M(r 



4.5 



that the normalizations 



Recall from ; 

of the FIRES cluster and CDFS structures are calculated as- 
suming a comoving line-of-sight extent of 10 Mpc. We also 
show the local m ass function of star- forming galaxies, from 
Peng et al. (20101, renormalized to minimize the x 2 value of 



the dense and underdense samples over the range for which 
the data are complete. This shows that both high density 
environments are a factor ~ 45 times denser than the field 
sample. This depends on our assumption about the line-of- 
sight extent. A firm lower limit to the overdensity is a factor 
7, obtained by assuming the cosmological volume between 
the redshift limits used to define each sample. We note that 
the shape of the single Schechter function, with M* and a 
fixed to their local values, is a good fit to both the field 
and overdense samples. There is no evidence for it to vary 
with environment, and the reduced x 2 is near unity for both 
samples. 



5.3 Specific star formation rate 

The sSFR of galaxies in our sample is shown as a function 
of their stellar mass in Figure |17| The small filled circles 
represent the individual emission line galaxies. The large 
symbols represent the binned mean sSFR of the present, 
combined sample (solid red circles), compared with ROLES 
at z = 1 (solid black circles), and local star-forming SDSS 
(blue diamonds) datasets. 

We find a distinct anti-correlation between sSFR and 
M*, with a power-law slope of /3 ~ —0.2 over the mass range 
10 9 < M*/M < 10 10 . This is similar to our finding at 
z — 1 with ROLES, and steeper than the anti-correlation 



at z — which already poses a challenge to models (Bower 



et al. 2012 Weinmann et al.||2012[ ). At all stellar masses, 
the sSFR at z = 0.75 and z = 1 is significantly larger than 
locally. Moreover, the relation significantly departs from a 
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Figure 17. The sSFR of the star-forming galaxies in our sample 
are shown as a function of stellar mass. Small points represent 
individual galaxies, while the filled circles with error bars are 
the mean values, in equal mass bins. Comparable mean values 
are shown at z = 1 from ROLES (solid black circles), and at 
< z < 0.1 from Stripe82 ( |Gilbank et al.||2011|>, w ithout error 
bars for clarity. We overplot the |Gilbank et al.| ( |2011| models (see 
text for description) at three redshifts, as indicated in the legend. 
Also shown is the mass limit of this survey (thick grey vertical 
dashed line), and the SFR limit (dashed, light blue line). 
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Figure 18. As Figure[l7] but where our data are now divided into 
high density structures, the MS1054 cluster (open green squares), 
and field environments. We show only the Gilba nk" et al.| (|2011) 
model corresponding to z = 0.75, as the smooth, red curve. 



simple power-law, with an upturn observed at low stellar 
mass (see also [Brinchmann et al.||2004 Elbaz et al.| |20Q7 



Popesso et al.|2011 ). To interpret this, we turn to the models 
used by Gilbank et al. (2011), based on the staged galaxy 
formation models of Noeske et al.l p007 ) . These are shown 



as the smooth curves on Figure |17| In this model, galaxies 
are parametrised by an exponentially declining SFR, with 
formation redshift and SFR timescales both a function of 
stellar mass. This simple description provides a reasonable 
match to the observations at all three epochs shown, in- 
cluding the increase in sSFR observed at the lowest stellar 
masses in both the present data and ROLES. 

Finally, in Figure [18] we show the sSFR-M* relation in 
different environments. Both high and typical density popu- 



lations show a decreasing sSFR with increasing stellar mass. 
The shape and normalization of the relation in all environ- 
ments are consistent with one another, and with the model 



of Gilbank et al. (2011), over the entire stellar mass range. 



6 DISCUSSION 

It has consistently been shown that the main influence of 
environment is on the fraction of star-forming galaxies of a 
given stellar mass (e.g. Bal dry et al.|2QQ6j |Wijesinghe et al.| 
2012 ). Amongst the star-forming population itself, any resid- 
ual environment dependence is weak at best (e.g. Peng et al. 
2010). However, a transformation from active to passive can- 



not be instantaneous and must therefore be reflected in an 
environmental dependence of the sSFR distribution. Results 
here are more controversial and depend on how the analy- 
sis was done. At low redshifts, z < 0.5, dense environments 
consistently show either no change in the sSFR relation (e.g. 



Peng et al. 2010) or a mildly reduced sSFR (e.g. Vulcani 
et al. 2010). At z ~ 1 and above, there is some evidence 
that sSFR is higher amongst star forming galaxies in some 
dense environments, but not all ( |Sobral et al.poTT ). While 
we will not resolve this issue here, our data extends the dis- 
cussion to low-mass galaxies at an intermediate redshift of 
z ~ 0.75. 

In Figure [19] we re-create the density segregated, 
sSFR-M* plot from |Li et aLl(2011| their Figure 6). In this 
plot we show the density-dependent sSFR— M* relation for 
each of the z = 0.1 (SDSS), z = 0.75 (present study) and 
z — 1 (ROLES) epochs. The SDSS and ROLES samples were 
segregated according t o a loca l densi ty parameter, p 5 , which 
is defined in detaip]in |Li et al. ( 2011 ). In the Figure, we show 
only the relations for the average environment, and the most 
overdense subsample, for best comparison with the present 
data. For the ROLES sample at z = l, lLi et al.|fl201lD find 
the local density ranges from 0.04 Mpc to 0.6 Mpc - , with 
an average of ~ 0.1 Mpc 3 . Thus their overdense regions are 
a factor ~ 6 denser than the average. We estimated the 
volume density of our subsamples in § |4.3| assuming a 10 
Mpc line-of-sight extent for our high density structures this 
is ~ 0.9 Mpc -3 , compared with the field value of ~ 0.02 
Mpc -3 . While these numbers cannot be compared directly 
with p5, they are not too dissimilar in practice. The small 
fields mean that the nearest-neighbour approach taken by 



Li et al. (2011) is almost exclusively identifying structure in 
redshift space, as we do here. The main difference is that Li 
et al. use massive galaxies as tracers, while we use low- mass 
emission line galaxies. While the physical scales over which 
the density is estimated will not be identical, we may expect 
that the relative density between structures and the field 
can be fairly compared between the two analyses, to within 
a factor of a few. Thus, the high-density environments of the 



4 Briefly, p§ indicates the redshift-completeness weighted num- 
ber density of star- forming galaxies limited to Mk ab ^ —21.0, 
found within a 'nearest-neighbour' volume defined by the five 
closest galaxies to the current galaxy being evaluated. The vol- 
ume is defined by the maximum projected radius of the set of five 
nearest-neighbour galaxies, and the difference in co-moving dis- 
tances set by the closest and farthest nearest-neighbour galaxies 
in redshift space. 
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Figure 19. The mean correlation between sSFR and stellar mass 
is shown for our sample, as the points with error bars. The pur- 
ple diamonds represent the result in the "field", while the red 
circles represent the high-density environments (MS 1054 cluster 
and CDFS structures), that are a factor ~ 45 overdense. This is 
compared with similar relations at z = 1 (ROLES) and z = 
(SDSS), taken from the analysis by |Li et al.| ( |2011| ) . In this case, 
the red line only represents an overdensity of ~ 6 relative to the 
average environment, represented by the purple line. The thick 
grey vertical dashed line highlights the mass limit of the survey. 



present study, overdense by a factor ~ 45, are likely to be 



significantly denser than those of Li et al. (2011 ). In fact, 45 
is probably an underestimate, as at least the MS 1054 virial 
diameter is likely much less than 10 Mpc. 

Considering first the average, "field" environments, Fig- 
ure [19] shows smooth evolution of the sSFR — M* relation 
from z = to z = 1, with little or no significant change in 
slope, but a normalization that increases approximately as 



In contrast, the low-mass end slope of the rela- 



tion in high-density environments shows mild evolution. It 
is natter than the average relation locally, and steeper than 
the average relation at z = 1. Interestingly, our new data at 
an intermediate redshift z ~ 0.7 shows no difference at all 
between the two environments. Since our survey includes en- 
vironments with much higher densities, the actual evolution 
between z ~ 0.7 and z ~ 1.0 of comparably dense envi- 
ronments could be considerably stronger than shown here. 
Note that the analysis of |Li et aL (|2011| includes an under- 



dense environment, for which the contrast with their densest 
environments is considerably larger. 

The implication is that the sensitivity of low-mass, star 
forming galaxies to their environment has evolved signifi- 
cantly from z = 1 to the present day. Today, the average 
sSFR of such galaxies is slightly lower in high-density envi- 
ronments, while at z = 1 the average is slightly higher. This 
"reversal" of the SFR-environment relation has been noted 
by others (e.g. Elbaz et al.|2007 Li et al.|2011 ), and our new 
data at z — 0.7 appear to correspond to the "transition" 
epoch where the sSFR-mass relation shows no environmen- 
tal dependence. 



7 CONCLUSIONS 

For the first time, the faint [Oil] emission from low stel- 
lar mass galaxies (8.5 < log (M*/M@) < 9.5) has been 
spectroscopically measured for galaxies in the redshift range 



0.62 < z < 1.15. By targetting fields (CDFS and FIRES) 
with known overdensities, including the massive MS 1054-03 
cluster, we explore how star formation in these low-mass 
galaxies are affected by their environment. Our main con- 
clusions are as follows: 

• There is little, if any evolution in the galaxy stellar mass 
function of [Oil] luminous galaxies between z = and z ~ 1. 

• The trend of a decreasing specific star formation rate 
with increasing stellar mass has been confirmed down to 
unprecedented stellar masses. The normalization at z ~ 0.75 
is similar to our earlier results at z = 1, and significantly 
higher than at z = 0. 

• The average power- lase of the sSFR— M* relation is /3 ~ 
—0.2, with indication of a steeper relation at low masses. 
This is consistent with what we found at z = 1 with ROLES. 

• The star formation rate density shows little evolu- 
tion between z = 0.7 and z = 1, but is consistent with 
the (1 + z) 2,0 evolution expected from comparison with 
the SDSS. The SFRD evolution is consistent with a mass- 
independent evolution in normalization; that is, the char- 
acteristic mass of star-forming galaxies is independent of 
redshift. However we caution that systematic and statisti- 
cal uncertainties preclude us from establishing the SFRD to 
better than a factor of ~ 2 at any mass; thus there is room 
for relatively small, mass-dependent evolution. 

• Environment is found not to influence the sSFR — M* 
relationship at any stellar mass at the epoch studied here, 
z ~ 0.75. This suggests that the SFR of star- forming galax- 
ies is not enhanced or diminished by local density, and that 
the apparent reversal in correlation between environment 
and star formation rate occurs at higher redshift. 

Our results on the constancy of the SF-GSMF are con- 
sistent with many previous works 



(e.g. 



Pozzetti et al. 2010 



Gilb ank et al.poTT ), extending the results to lower masses 
at this intermediate redshift. In Li et al. (2011) we showed 



the emergence of a weak environmental dependence of the 
sSFR-mass relation on environment, such that the lowest 
mass galaxies in the densest environments show a signif- 
icant excess sSFR compared to their lower density coun- 
terparts. In this work, we show the absence of this trend 
~1 Gyr later. By the present day, another 5 Gyr later, the 
relation has reversed, and low-mass galaxies in dense envi- 
ronments have lower sSFR than the average. This smooth 
transition is consistent with e.g. IQuadri et al. (2012), who 



pointed out that it would be strange to see a sharp transi- 
tion given factors such as the apparently smooth growth of 
passive galaxies over cosmic time. We note, however, that 
this environmental-dependence of the star-forming popula- 
tion is very mild, and much weaker than the aforementioned 
evolution in the fraction of galaxies with no star formation 
whatsoever. 
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